i£> 

© 

©0 

in 


R-2209-ARPA 
March  1978 


ARPA  ORDER  NO.:  3486 
8G10  Tactical  Technology 


D D C 

ESIpnOilE 

APR  24  I97S 


EtTFE 

B 


Simple  Relations  for  the  Stability  of 
Heated  Laminar  Boundary  Layers  in  Water: 
Modified  Dunn-Lin  Method 

J.  Aroesty,  W.  S.  King  G.  M.  Harpole, 

W.  Matyskiela,  A.  R.  Wazzan,  C.  Gazley,  Jr. 


A Report  prepared  for 

DEFENSE  ADVANCED  RESEARCH  PROJECTS  AGENCY 


"‘WSTWBUflON  STATEMEWtX 

,'.y  

approved  far  public  mImmw 

* DlitilbutioB  Unlimited 


Rand 

SANTA  MONICA  CA  SOW 


The  research  described  in  this  report  was  sponsored  by  the  Defense  Advanced 
Research  Projects  Agency  under  Contract  No.  MDA903-78-C-0121 . 


Reports  of  The  Rand  Corporation  do  not  necessarily  reflect  the  opinions  or 
policies  of  the  sponsors  of  Rand  research. 


PubtUiied  by  The  Rand  Corporation 


Qsi  ^ ^ ARPA  QRDE £bI^  3486 
8G10  Tactical  Technology 


/R-2^9-ARPA 
m.  /:r7l  Manrfi»78 


IIj  A L 


D D C 

aarnoiiE 

APR  24  1978 

EGSEDirE 

B 


Simple  Relations  for  the  Stability  of 
Heated  Laminar  Boundary  Layers  in  Water: 
- — Modified  Dunn-Lin  Method  . i 


./Aroesty,  W.  S./Kin^  G M/^HarpoN 


A.R. 


A Report  prepared  for 

DEFENSE  ADVANCED  RESEARCH  PROJECTS  AGENCY 


Randi 

MONICA 


1% 


AmOVEO  FOi  FViilC  RELEA&£i  tM»ft(tUTiON  UMtUUf  EO 


-iil- 


PRKFACE 

Under  the  sponsorship  of  the  Tactical  Technology  Office  of  the 
Defense  Advanced  Research  Projects  Agency,  Rand  has  been  investigating 
Che  fluid  mechanics  and  hydrodynamics  of  low-drag  submersible  vehicles. 

Preliminary  studies  of  such  vehicles  often  involve  characterization 
of  surface  heating  effects  on  the  hydrodynamic  stability  of  lamitiar 
boundary  layers  in  water.  Elaborate  numerical  solutions  of  the  Orr- 
Sommerfcld  equations  are  usually  required  for  such  characterization, 
even  for  preliminary  conceptual  studies. 

An  analytic  method  first  developed  more  than  30  years  ago  is  re- 
viewed and  modified  to  provide  a simple  set  of  relations  for  determin- 
ing the  effect  of  heating  and  pressure  gradient  on  the  minimum  critical 
Reynolds  number  of  laminar  boundary  layers  in  water.  This  new  method 
does  not  obviate  the  need  for  numerical  solution  of  the  Orr-Soramerfeld 

9 

equation,  such  as  might  be  required  for  "e  " calculations.  However, 
it  does  provide  an  inexpensive  and  simple  way  to  estimate  the  effects 
of  temperature  and  pressure  gradient  on  the  minimum  critical 

Reynolds  number. 

This  report  should  be  useful  to  hydrodynamic ists.  designers  of 
submersibles,  and  others  engaged  in  the  application  of  fluid  mechanics 
to  the  improvement  of  underwater  vehicle  performance.  Related  Rand 
reports  include: 
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SUMMARY 

A simplified  method  is  described  for  calculating  the  effect  of 
surface  heating  on  the  hydrodynamic  stability  of  heated  laminar 
boundary  layers  in  water.  The  method  involves  modification  and  up- 
dating of  relations  first  developed  by  Lin  for  a constant-property 
fluid  and  later  by  Dunn  and  Lin  for  a compressible  one.  These  modi- 
fied Dunn-Lin  relations  for  the  minimum  critical  Reynolds  number  give 
results  in  substantial  agreement  with  numerical  integration  of  the 
Or r- Sommer f eld  equations.  The  method  is  then  employed  to  evaluate 
the  influence  of  ambient  temperature  level  and  wall  temperature  dis- 
tribution on  the  minimum  critical  Reynolds  number.  Preliminary 
studies  indicate  that  wall  temperature  distributions  can  be  found 
that  have  a strong  favorable  effect  on  the  stability  of  laminar 
boundary  layers  in  water,  even  for  adverse  pressure  gradient  flow. 
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I . INTRODUCTION 


Successful  laminar  flow  hydrodynamics  relies  on  the  manipulation 
of  velocity  profiles  in  the  wall  region  of  the  laminar  boundary  layer 
to  enhance  flow  stability.  Pressure  gradient,  body  shaping,  suction, 
and  surface  heat  transfer  are  classic  methods  for  accomplishing  this 
manipulation. 

Recent  theoretical  and  experimental  studies  suggest  that  surface 
heating  in  water  holds  promise  for  enhancing  boundary-layer  stability 
by  slowing  the  growth  of  Tollmeln-Schlichting  instabilities,  or  by 
increasing  the  surface  area  over  which  two-dimensional,  infinitesimal 
disturbances  are  damped.  In  the  absence  of  a comprehensive  theory  of 
boundary-layer  transition,  linear  stability  theory  currently  provides 
the  sole  analytic  guide  for  manipulating  mean  flow  velocity  profiles 
to  delay  transition. 

Even  for  preliminary  engineering  applications,  no  single  parameter 
describes  the  stability  of  a particular  velocity  profile.  However, 
the  minimum  critical  Reynolds  number  can  be  extremely  useful,  both  as 
a precise  measure  of  the  extent  over  which  all  two-dimensional,  infin- 
itesimal disturbances  are  damped  and  as  a simple  qualitative  surrogate 
for  the  stability  characteristics  of  a particular  velocity  profile  [1]. 

In  1946,  C.  C.  Lin  (2)  published  a simple  set  of  relations  for 
determining  the  minimum  critical  Reynolds  number  from  the  velocity  pro- 
file of  a constant  property  laminar  boundary  layer.  The  original  Lin 
relations  (2,3], 


v(c)(l  - 2X)  = .58  , 


(la) 


•k 

Recent  computations  by  Wazaan  and  Gazley  (1]  suggest  that  the 
minimum  critical  Reynolds  number  may  also  be  usieful  in  correlating 
transition  predictions  obtained  by  the  "e^”  method. 
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* 
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where 


(lb) 


(Ic) 


25  1^  (0) 

^^^6^rain-crit  4 

c 


(Id) 


= height  of  critical  layer,  and  c = wave  speed,  normalized  to  free- 
stream  velocity,  were  based  on  asymptotic  analysis  of  the  Orr-Sommerfeld 
equations  (valid  for  large  values  of  cxRe^,  small  values  of  the  wave- 
length, a,  and  small  values  of  c)  and  a numerical  factor  derived  from 
a more  complete  calculational  method.  These  relations  were  widely 
used  (31  until  the  advent  of  computer-based  schemes  for  the  solution 
of  the  Orr-Sommerfeld  equations.  In  1946,  Lees  and  Lin  (4)  extended 
the  original  Lin  relations  to  the  compressible  flow  of  air,  aivl  later, 
Dunn  and  Lin  (5)  and  Mack  [61  presented  further  compressible  extensions 
of  the  original  Lin  analysis  of  the  neutral  curve  separating  stable  and 
unstable  regions.  However,  special  difficulties  associated  with  com- 
pressibility diminished  the  role  of  this  analytic  approach  in  high- 
speed dynamics,  and  computer-based  numerical  solutions  have  since  pre- 
empted stability  analysis  for  air  flows. 

For  the  case  of  practical  heated  water  boundary  layers,  the  situa- 
tion is  simpler.  Water  density  is  nearly  constant,  temperature  and 
viscosity  fluctuations  have  little  effect  on  stability,  and  the  primary 
departure  from  the  constant  property  ineompressible  flow  originally 
considered  by  Lin  is  the  variation  of  mean  flow  viseosity  with 
temperature. 

On  the  basis  of  numerleal  testing  Uaszan  (7,B]  suggested  that 
viscosity  and  temperature  fluctuations  are  negligible  and  tliat  the 
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appropriace  Orr-Sommcrfeld  equation  for  heated  water  boundary  lave. s 
Is  still  fourth  order  but  includes  derivatives  of  mean  flow  viscosity 
with  normal  distance.  The  equation  is 


(U  - cXr  - Ci(^) 


- 2a^0"  + 0*0 


- * (4)  t 

+ 2u' (<!>'•'  - aV')  + y'’(^"  + “ 0 , 


(2) 


where  primes  refer  to  d/dn;  H =*  Y/6;  Rg  = U6/v^;  the  disturbance  stream 

^ Cctx  “ c ti ) 

function  ilj(w,y,t)  = <Xn)e  ; and  o and  c are  real  along  the 

neutral  curve.  Both  velocity  and  viscosity  are  normalized  by  their 
values  as  n ’’  ‘"o  It  was  also  demonstrated  that  the  omission  of  the 
terms  u'  and  y"  had  a small  effect  (~  3%)  on  the  numerical  calculation 
of  minimum  critical  Reynolds  number.  This  is  not  surprising  for  the 
extremely  stable  flows  sometimes  encountered  with  heated  water  boundary 
layers,  because  the  derivatives  of  are  large  in  the  critical  region, 
and  the  terms  y*  and  y"  are  essentially  of  unit  magnitude  or  less. 

In  this  brief  report,  we  show  how  the  Lin  relations  and  the  Dunn- 
Lin  theory  may  be  used  after  slight  modification  to  estimate  the  min" 
imum  critical  Reynolds  number  for  heated  water  boundary  lay^Jts.  The 
appropriate  approximations  and  dimensional  factors  in  these  new  rela- 
tions are  based  on  recent  calculations  of  the  Orr-Sommerfeld  equations 
for  constant  property  flow.  The  results  of  our  modified  bunn-Lin  theory 
are  then  compared  with  values  obtained  from  numerical  studies  of  the 
Orr-Sommerfeld  equation  for  heated  water  houndary  layers.  On  ? the 
accuracy  of  these  new  relations  has  been  established,  they  are  used  to 
indicate  the  effect  of  changing  temperature  levels  on  the  minimua  crit- 
ical Reynolds  number.  The  surprisingly  powerful  effect  of  te«p«’’‘ature 
variation  in  the  maintenance  of  a stable  flow  is  illustrated  by  the 
case  of  an  unfavorable  pressure  gradient,  because  of  r'-'ir  accuracy 
and  simplicity,  these  relations  can  be  useful  in  preliminary  feasi- 
bility and  optimisation  studies  of  laminar  flow  hydrodynamics. 


It.  REVIEW  AND  ANALYSIS 


?;>r  jnv  .>01  ar'rc,  we  briefly  review  the  Duim-Lin  method  for  solving 
i^q.  (2)  along  Che  ne<  crnl  curve  (but  only  chose  aspects  of  the  method 
chat  are  InstrumenCai  in  this  modification;  complete  elaboration  is 
found  in  Uefs.  2,  3,  or  6). 

Because  aR^  is  e»pecCed  to  be  large,  classic  asymptotic  methods 

should  be  sufficiently  accurate  for  engineering  estimates.  I''ir.st,  the 

inviscid  Urr-Sommcrfeld  equation  obtained  after  neglecting  terms  of 

0(l/aRo)  is  considered: 

0 

(U  - cXr  - a(J))  - U"(D  - 0 , (3) 

Following  Lin  (2],  the  solution  of  this  equation  that  satisfies  the 
outer  boundary  conditions  0(®)  » 0'(“*)  » 0 is  called  (^(n)  and  is  cal- 
culated by  a power  series  in 

The  viscous  solution  is  more  complicated.  The  modified  Tollaien 
variation,  C,  used  for  obtaining  the  dominant  viscous  solutions  of 
Eq.  (2),  accounts  for  both  departures  from  linearity  of  the  velocity 
profile  and  variable  viscosity.  It  is  defined  by 


Tlte  leading  terms  of  £q.  (2)  then  become 

- iC  f.AO  0 . (%) 

where  F(C)  is  ttw  doeinant  viscous  solution  of  the  Orr-Soeiaerfeld 
equation.  There  are  fou;  independettt  solutions  of  this  equation,  two 
are  essentially  inviscid  solutions  and  are  already  included  in  the 
inviscid  solution  $Crt)>  and  one  et  tite  reoaining  solutions  is  unbounded 
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asj  ^ ■*  ».  Thus,  there  ts  only  one  acceptable  viscotis  solution.  This 

solution  8<ust  be  multiplied  by  a factor  to  give  the  correct  asyraptotir 

forw  as  aR(5  ■»  *.  This  factor  introduces  an  error  not  larger  than 
-1/3 

O(aR^)  ' in  Eq.  (b) . 

The  appropriate  viscous  solution  is  ^ , (n) > where 


4>vi^(n)  - - i(aR^) 


-1/2 


a{ti) 


(6a) 


•/ou 


aiul 


A(n) 


^,-(n)/3 


(6b) 


Thus  the  required  approximate  solution  is  the  linear  combination 

a(|i(n)  * 

The  no-alip  conditions  at  the  wall,  i^(0)  « (^*(0)  » 0,  must  still 
be  satisfied.  This  leads  to  the  relations 


mo)  4.  » 0 . 

(?a) 

a4>(a)  4.  bf  . (0)  » 0 . 
y\m 

which  then  result  in  the  eigenvalue  relation 

(7b) 

$((» 

m 

m)  r.^(a)  ‘ 

where  i.  » 
vis 


4 aad 

4 ^*Cn). 
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The  ratio  of  the  viscous  solutions  can  then  be  written  as 


'vis»>  i “ 1^0’'“ 


c - U 


dn  F(Z)  , 


where 


(otRg) 


1/3  3 


and  F(Z),  the  Tietjen's  function,  is  defined  by 


F(Z)  = 


/ 

•'co  •'00 


H<J>  [f  au^«] 


- Z / hJJ*  [I  d? 


If  ^ is  defined  by 


Eq.  (9a)  becomes 


W!^ 


5f(by  • 

dn 


and  the  eigenvalue  relation  (£q.  (8))  becomes 


Max. . P(Z\ 

iMo)  u*(0)  * 
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It  is  traditional  to  write  this  equation  in  the  form 

E(a,c)  = F(Z)  , (13a) 


where 


E(a,c) 


<t»(0)  U*(0) 

Hj;(0)  c(l  + X)  * 


(13b) 


The  quantity  X is  usually  small.  If  both  velocity  gradient  and 
viscosity  are  constant,  the  X is  identically  zero.  In  other  more 
realistic  cases,  departures  from  this  constancy  in  the  region  between 
wall  and  critical  layer  are  generally  small  enough  to  result  in  jX]  <<  1. 
However,  its  selective  inclusion  in  numerical  calculation  results  in 
inore  accurate  estimates. 

Equation  (13a)  can  be  recovered  from  the  original  Dunn-Lin  formu- 
lation by  setting  '•*  0 and  = 1 and  permittijig  y to  vary  with  r|. 

It  is  also  customary  to  introduce  the  modified  Tietjens  function  de- 
fined by  3(Z)  = 1/(1  - F(Z)1,  and  to  rewrite  E(a,c)  in  terms  of  auxiliary 
real  valued  functions  u(a,c)  and  v(a,c)  defined  by 

(1  X)(u  iv)  . .. 

^-Vc) 

The  .secular  relation  (Eq.  (13a))  then  becomes 

’ 3(Z)  » 3^(2)  + i3j(Z)  , (13) 


(1  4 X)  (u  iv) 
I + X(u  iv) 


where  3(Z)  is  decomposed  into  real  and  imaginary  components. 

Lin  (2)  showed  ttiat  the  inviscid  solution  could  be  written  as 


u ♦ iv  » 1 + u'<0)c 


(I  - e)"^  ♦ (1  - c)"“  K,®*  ♦ KjCt’  ♦ . . . 

a *-  Ci  - c)’^ti^a^  ♦ - c)"^  ^ . . . 


» 


m) 
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c 

i 


where  K^(c)  = | (U  - c)  ^ dy  , 


' I 

/ 


c y 

-2  . /•  .2 


K2(c)  =1  (U  - c)  dy  / (U  - c)  dy  , 


I 


' I 

/ 


c y y 

-2  . /•  x2  . f .-2 


^3(0)  = I (U  - c)  dy  I (U  - c)  dy  I (U  - c)  dy  , 


L 


(17) 


t 

Hj^(c)  =»  I 
J (\ 


(U  - c)^  dy  , 


/■ 


H2(c)  - / (U  - c)^  dy  /"  (U  - c)^  dy 


> 

i 


n.. 


( 


j( 


HjCc)  » I (U  - c)^  dy  ^ (U  - c)  ^ dy  ^ (U  - c)^  dy  . 

Jo 


By  direct  evaluation  of  Eqs.  <16)  and  (1?)  tor  i»talX  a and  c,  Lin 
found  that  Eq.  (lb). 


wU'(0)U(n  )U*’(n  ) 

V V 


is  a suitable  approx iaat ion  for  v. 

Lin  also  detensined  titat  the  Ainiautt  value  of  Ke^  is  opproxiaated 
when  a^(S)  has  its  aaxitaua  value,  correspv^^nding  to 
• .SB,  * l^i^^B.  and  7.^  > 3.21.  He  then  propoeed  on  Iteration 


t* 
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procedure,  valid  for  small  values  of  A,  for  solving  Eq,  (15).  The 
procedure  is  based  on  neglecting  terms  of  0(X)  in  evaluating  U,  but 
retaining  such  terms  in  the  evaluation  of  v.  From  Eq.  (15): 


u E!  u (Z,,) 

r C 


(18a) 


|n-x[i  - 23^(Z„)]}=  a^(y 


(18b) 


Equation  (la),  Lin’s  original  relation  for  determining  c,  thus 
fv-x'iows  from  Eq.  (18b).  The  relation  between  oRe^  and  Z follows  from 
Eqs.  (9b)  and  (10): 


'“““{'critical 


a)'  u'(o)^  u( 

3 2 

c-*(l  9-  X)^ 


Originally,  Lin  estimated  cs  from  a crude  numerical  solution  of  Eq. 
(16)  for  the  hlasiusi  profile, 

a Sf  li’  (0)c  . ' 

He  thfc.1  approximated  Eq.  (19)  by 


(It.,) 

0 critical 


after  neglecting  X and  adjusting  his  constants  to  agree  with  results 
obtained  fr<»B  his  own  sore  detailed  calculations.  In  the  saae  spirit 
as  Lin,  we  also  assume  that  ei  « U*(0)e,  neglect  terms  of  9(X)  in  Eq. 
(19).  and  deduce 


(S,  > 28  ()•«))  U(0) 

It  U-al  “ ^,4  , 
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where  the  numerical  constant  is  now  chosen  to  agree  with  recent  accu- 
rate numerical  calculations  for  the  constant  property  Blasius  boundary 
layer. 

Because  1 is  a small  quantity,  we  bypass  the  required  numerical 
quadrature  between  wall  and  critical  layer,  and  with  sufficient  accu- 
racy estimate  it  by 

This  formula  results  from  the  assumption  of  a linear  viscosity  profile 
and  parabolic  velocity  profile  between  wall  and  critical  layer  and  the 
evaluation  of  the  dominant  term  in  the  integral  in  £q.  (10).  The  re- 
quired modifications  of  the  Dunn-Lin  theory  are  then  Eqs.  (la),  (lb), 

(10),  and  (22)  for  determining  n and  c;  and  £q.  (21)  for  the  evalua- 

c 

tion  of  tl»e  minimum  critical  Reynolds  number, 


-11 


III.  RESULTS 


FIGURE  1 

A comparison  between  values  of  calculated  by  the 

present  approximate  method  and  those  calculated  by  numerical  solution 
of  the  Orr-Sommerfcld  equation  is  shown  in  Fig.  1.  The  lines  are  based 
on  published  results  of  Wazzan  and  co-workers  [7-‘101  for  Falkner-Skan 
similarity  flows  at  various  surface  overheating.  The  symbols  corre- 
spond to  the  current  results.  Although  no  single  profile  shape  factor 
can  uniquely  correlate  when  both  pressure  gradient  and 

surface  heating  are  present,  U = 6*/Q  remains  a convenient  parameter 
for  presenting  these  results. 

The  open  circles  are  computed  by  our  approximate  formulas  for 
wedge  flows  at  zero  AT.  The  agreement  between  exact  calculation  and 
the  proposed  method  is  remarkable  in  Che  practical  range  of  U between 
2.2  and  1. 

The  asterisks  correspond  to  the  case  of  S = 0 (flat  plate)  and 
various  surface  overheating.  Note  that  tl\e  agreement  between  exact 
and  approximate  methods  is  again  satisfactory  in  Cite  region  between 
H e 2.59  and  H « 2.2,  where  a stability  reversal  occurs  at  AT  Increases 
past  75*F.  The  method  agrees  with  exact  calculations  in  predicting  tlie 
existence  and  location  of  the  maximum  attainable  critical  Reynolds 
number,  both  for  this  case  and  for  the  case  of  6 * -.1988  (closed 
circles) . 

For  8 = 0,  the  method  appears  to  overpredict 
the  region  of  large  surface  overheating  corresponding  to  values  of  H 
less  than  2.1. 

Tlie  case  of  8 = l.O  (crosses),  corresponding  to  an  extremely 

favorable  pressure  gradietit,  is  more  problematic.  F.xact  calculations 

. 4 

indicate  tlrat  (Re.^)  , exhibits  a local  maximum  at  2.2  s lO  and 

0"  critical 

U = 2.05  and  then  increases  again  with  further  surface  overheating, 
ttie  approximate  method,  however,  exhibits  an  absolute  maximum  at 
U = 2.05,  corresponding  to  AT  « 40*F.  Tte  agreement  between  ttie  exact 
and  approximate  values  is  good  up  to  this  value  of  surface  overheating. 


Local  max  at  AT~40® 


work 


“ Uothcrmol  wedge  flows^ 
Ref.  9 

Wedge  flows  with  hoot  transfer 
in  woter,  a 

Q w -0.1968  Ref.  7 

^ * 0 Ref.  10 

Q - 1.0  Ref.  8 
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The  modified  Dunn-Lin  method  can  be  used  with  reasonable  con- 
fidence when  moderate  surface  overheating,  of  the  type  achievable 
in  practice,  is  considered. 

FIGURE  2 

Figure  2 shows  the  varlstlon  of  critical  surface  over- 

heat for  a series  of  favorable  and  unfavorable  pressure  gradients  of 
the  Falkner-Skan  type.  The  symbols  in  this  figure  represent  results 
obtained  by  numerical  integration  of  the  Orr-Sommerfeld  equations, 
and  the  lines  refer  to  the  approximate  method. 

The  agreement  between  exact  and  approximate  computations  is  again 
satisfactory  in  the  practical  range  of  "Bs"  between  +.10  and  -.10. 
These  results  support  our  previous  observation  about  the  method  being 
accurate  for  practical  levels  of  pressure  gradient  and  surface  over- 
heating. 


FIGURE  3 

In  Fig.  3,  the  modified  Dunn-Lin  method  is  used  to  calculate  the 
dependence  of  aohient  water  temperature  for  several 

surface  overheats  and  Bs  of  .10,  0,  and  -.10,  The  variation  of  water 
viscosity  with  temperature  is  the  primary  mechanism  for  the  surpris- 
ingly large  effect  of  ambient  water  temperature  on 
corresponding  to  a three-fold  increase  when  ambient  temperature  is 
decreased  from  7S*F  to  4S*F. 


FIGURE  4 

All  of  the  previous  results  are  for  constant  w.%11  temperature. 
Figure  4 compares  results  obtained  by  the  present  method  with  exact 
computations  for  variable  wall  temperature  (AT  X)  and  .m  adverse 
pressure  gradient  corresponding  to  B • -.10.  The  ability  of  the 
approximate  method  to  predict  the  qualitative  features  of  the 
^^§*^erttical  again  demonstrated.  Sots  that 

the  variable  wall  temperature  increases  stability  over  the  constant 
wall  temperature  case  wfien  the  local  temperature  difference  is  greater 
than  4S*F.  The  surprisingly  poserful  effect  of  wall  temperature 


Modified  Dunn -Lin 
approximation 

/ 

/f 
// 
/ k 

/\ 

^ Exact  computotions 

h 
/ / 

^AT-x 

f 

0=  -0.10 

/ 1 

1/ 

0 

>*AT  = Constont  (exact 
/ computations  and 

present  method  give 
similar  results) 

Locol  ^T«  *F 

FIq.  4'~TK«  tfftct  ef  voffobt*  lurfow  otf«fK«o>  on 
nouteol  itobility:  Modifiod  Oonn^'Lin  vt.  oiUKt 
cOHiy>ulolions  • 6^*) 

(Not*  tht  incrooM  in  ifobility  ovor  tho  coratont 
vmll  loMotfolufo  COM  > 4S*F) 
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varlation  in  Increasing  the  extent  of  stable  laminar  flow  is  due  pri- 
marily to  the  upstream  history  of  the  thermal  boundary  layer  and  its 
ability  to  appropriately  deform  viscosity  and  velocity  profiles  in 
the  region  between  critical  layer  and  wall. 

CONCLUSION 

W«  have  developed  a modification  of  Lin's  method  that  can  be 
used  to  estimate  the  minimum  critical  Reynolds  number  of  a heated 
laminar  boundary  layer.  This  method  is  both  simple  and  surprisingly 
accurate.  As  in  Lin's  original  work,  the  method  is  most  applicable 
when  c 0,  a 0,  and  ■*  “.  The  internal  details  of  the  method 

do  not  agree  entirely  with  exact  analyses  in  other  ranges  of  these 
parameters,  but  the  few  required  constants  have  been  chosen  so  that 
predictions  of  minimum  critical  Reynolds  number  are  consistent  with 
results  of  more  elaborate  numerical  integration  of  the  Orr-Sommerfeld 
equatious. 

This  method,  which  we  have  labeled  the  "modified  Dunn-Lin  method," 
can  be  used  with  confidence  In  engineering  studies  of  laminar  boundary 
layer  control  in  water.  Preliminary  studies,  based  on  this  method, 
suggest  that  an  appropriate  suriace  temperature  distribution  can  promote 
stability,  even  for  flows  experiencing  adverse  pressure  gradients. 


19 


REFERENCgS 


1.  A.  K-.  Wazzcin  and  C.  Gazlcy,  Jr.,  "The  Combined  Effects  of  Pressure 

Gradient  and  Heating  on  the  Stability  and  Transition  of  Water 
I^oundary  Layers,"  Pvooeedinga  2>id  Intepnatioml  Confevenae  on 
Di'ag  Reduction,  London,  19??;  also  T!ie  Rand  Corporation, 
R-2175-ARPA. 

2.  C.  C.  Lin,  "On  the  Stability  of  IVo-Dimensional  Parallel  Flows," 

Quart.  Applied  Math.,  No.  4,  Vol.  Ill,  Part  HI,  January  1946, 
pp.  277-101. 

3.  C.  C.  Lin,  The  Theory  of  Uydrodynaitia  Stability,  Cambridge  Uni- 

versity Press,  London,  1966. 

4.  L.  Lees  and  C.  C.  Lin,  Investigation  of  the  Stability  of  the 

Laninav  boundary  Layer  in  a Compveeoible  Fluid,  NACA  TN  No.  1115, 
Washington,  D.C.,  September  1946. 

3.  U.  W.  Dunn  and  C.  C.  Lin,  " On  the  Stability  of  tlie  l.amlnar 

Boundary  Layer  in  a Compressible  Fluid,"  J.  Aeroap.  Sot.,  Vol. 

22.  July  1955,  pp.  455-477. 

6.  L.  M.  Mack,  "Numerical  Calculation  of  the  Stability  of  the  Com- 

pressible Laminar  Boundary  Layer,"  Jet  Propulsion  Laboratory, 
Report  No.  20-122,  May  I,  I960. 

7.  A.  R.  Wazzan,  T.  T.  Okamura,  and  A.  M.  0.  Smith,  "The  Stability 

and  Transition  ot  Heated  and  Cooled  Incompressible  Laminar 
Boundary  Layers,"  i’wo.  International  Meat  Tvanaf&r  Conference, 
FarCa-Veraaillea,  Vol.  II,  Paper  FC-1.4,  1970. 

8.  A.  R.  Wazzan,  C.  Keltner,  T.  T.  Okamura,  and  A.  M.  0.  Smith, 

"Spatial  Stability  of  Stagnation  Water  Boundary  Layer  with 
Heat  Transfer,"  Fhya.  FI.,  Vol.  15.  1972,  pp.  2U4-2U6. 

9.  A.  8.  Wazzan,  T.  T.  Okamura,  and  A.  M.  0.  Smith,  "Spa'tal  and 

Temporal  Stability  Cliarts  for  the  Falkner-Skan  Boundary  Layer 
Profiles,"  Douglas  Aircraft  Company,  Report  No.  DAC-67086, 
September  1968. 

10.  A.  R.  Wazzan,  T.  T.  Okamura,  and  A.  M.  0.  Smith,  "The  Stability 
of  Water  Flow  over  Heated  and  Cooled  Flat  Plates,"  S.  Heat 
Trana.  AIM',  Series  C.  Vol.  90,  1968,  pp.  109-U4. 


s^f-ufo^v” rr»s»,'^'c*^iow,  ar  t.?-^  ix^r  'wxw  f»»i«  !'«»•»»«») 


>*  • 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
DEFORE  COMPLETING  FORM 


i COVT  ACCKUlOW  NOJ  1 AtCl^tAMr'}  CATakOfi  NUWACM 


» riTkf  imnii 


Sinrple  Relations  for  tlia  Stability  of  Heated 
Laminar  Bomdary  Layers  in  Water:  Modified  ^ 
Dunn-Lin  Method 


i TvFC  Of  liei>ONT  • PCKtoe  COVtACO 

Interim 


*.  ■ACMfOAMIMa  OMO.  MCAOMT  MUMKCN 


f *uTMOAf«) 


J.  Aroesty,  W.  S.  King,  G.  M,  Harpole, 

W.  Matyskiela,  A.  R.  Wazzan,  C.  Gazley,  Jr. 


k.  coMTki.ev  om  auMr  numaca^*; 


MnA903-78-C-0121 


) I COMTMOkLIMC  Okkicc  MAMC  kMO  AOOMCtS 


I PeAroAMiN&  OAfikNi:  ATiOn  NAMC  AMO  AOOAt&S  FHOGAAM  tu  AM£|i,T.  AAOjEC*.  TAiA 

_ ^ _ AACA  A *OAA  UMl^  MUMiC  A* 

Ihe  Rand  Corporation 

1700  M^dn  Street  y 

Santa  Monica,  California  90406 


»2.  AEAOAt  OATt 

March  1978 


I*.  NUMACA  Ok  AAOCt 

19 


A MOMITOAihO  AOCnCV  mAMC  A AOOACSV<kk<MAMA«  Mam  CsMmMiiA  OMcaJ  I>-  •CCUA|T«  CL*lt  r«/ >At* 

unclassified 


i>A  OKCLAIftklCATlOM  OOAMGAAO«mO 
tCMCOUkC 


t«  Ott'AIAuTlOM  ITATCmAMV  rWMMc  AaapM; 


i^proiAgd  for  Public  Release;  Distribution  unlimited 


tT  ITATCmCMT  (Qftk0  !•  49*  tf 


li  uA  •x^  a 40«AC/I|r  aMPlg^J 


Laminar  Bounds 
Nydbe^'namics 
Fluid  feidhanies* 
Subfi^ines 


underwater  Veluclee 


M AHaT  A4C  * fgoPWAtf  Hi«Af  A»  AAAk  A I 


* 401^ 


see  re^^erse  side 


unclassified 


Mevatw  CLAMtnCAVMjM  »'%««  ^AWI  tfii*  f «A<A»- 


SCCuRjty  cl  ASf.rtCATTOK  This  l»AOCrwhwi  Dsumnfmm) 


A sijiplified  nethod  for  calculating  the  effect  of 
surfaoE  heating  on  the  hydroc^amic  stability  of 
heated  laminar  boundary  layers  in  water.  Ihe 
nethod  involves  modificaticn  and  ipdating  of  relations 
first  developed  by  Lin  for  a oaistant  property 
fluid  and  later  by  Dunn  and  Iln  for  a oorpressible 
one.  'These  modified  Dunn-Lin  relations  for  the 
minimun  critical  Reynolds  nutrber  give  results  in 
substantial  agreement  with  numerical  integration 
of  the  Orr-Soimerfeld  equations.  'The  powerful 
effects  of  surface  benperature  distribution  on  the 
stability  of  laminar  boundary  layers  in  water  is 
illustrated,  even  for  adverse  pressure  gradient. 

Ref.  (author) 


Tfticlassified 

MCUHITV  CLAISiriCATlON  Of  THIS  PAOtflHun  En(«r»<<) 


